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ABSTRACT 


vTwo  methods  are  described  for  measuring  interferometric  fringe  contrast 
with  sensing  arrays.  The  absence  of  air  turbulence  allows  contrast  to 
be  determined  quickly  by  measuring  the  intensity  of  the  fringe  pattern 
at  phases  90°  apart.  The  presence  of  turbulence  requires  that  the  peak- 
to-peak  fringe  intensity  be  found  by  sampling  the  fringes  as  they  are 
scanned  across  the  array.  An  analysis  is  presented  showing  the  relation 
between  maximum  error  in  fringe  contrast  measurement  and  the  number  of 
samples  taken.  Experimental  results  are  presented  which  show  the  effect 


of  turbulence  on  fringe  contrast  measurement. 
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COHERENT  IMAGE  FRINGE  CONTRAST  MEASUREMENT 
WITH  SENSING  ARRAYS 
William  W.  Seemuller* 


INTRODUCTION 

The  technology  of  solid-state, light-sensitive, arrays  is  sufficiently  ad¬ 
vanced  at  this  time  to  make  these  devices  useful  in  a  variety  of  optical 
measurement  applications.  Initially,  arrays  were  useful  only  for  low 
resolution  imaging  applications  because  of  small  element  count  and  poor 
dynamic  range.  Presently,  arrays  having  resolution  equivalent  to  a  videcon 
with  a  dynamic  range  of  103  are  available.  This  performance  opens 
up  applications  for  arrays  in  the  areas  of  high  resolution  image  digi¬ 
tization,  densitometry,  and  image  mensuration.  The  chief  advantage  of 
the  array  as  an  optical  transducer  for  these  applications  (as  opposed  to 
a  videcon  or  image  dissector)  is  the  fixed  and  accurate  element  geometry 
which  insures  a  fixed  and  known  spatial  relation  between  pixels  within 
the  array  area.  Because  the  array  has  many  elements  and  because  it  is 
scanned  electronically,  it  is  also  a  very  rapid  device  and  eliminates  the 
need  for  mechanical  scanning. 


♦William  W.  Seemuller  is  employed  by  the  Research  Institute,  U.S.  Army 
Engineer  Topographic  Laboratories,  Fort  Bel  voir,  Virginia  22060 


This  paper  treats  the  problem  of  using  arrays  to  measure  the  contrast 
of  fringes  produced  by  an  interferometer,  when  the  fringe  pattern  Is  varied 
in  time  across  the  array,  by  varying  the  path  length  of  one  arm  of  the  inter¬ 
ferometer.  Although  this  requirement  arose  because  of  the  Heterodyne  Optical 
Correlator  (HOC),  the  work  described  is  applicable  to  any  situation  where  the 
accurate  measurement  of  fringe  position  and  contrast  is  required,  such  as  in 
optical  component  testing  [1].  The  HOC  was  developed  to  provide  a  rapid  means 
for  generating  elevation  contours  from  stereo, aerial -photograph  pairs  for 
map  production.  The  HOC  has  been  described  elsewhere  [2]  and  will  not  be 
discussed  in  detail  here.  Basically,  the  HOC  is  an  interferometer  in  which 
one  member  of  a  stereophotograph  pair  is  inserted  in  each  arm.  The  two  photo¬ 
graphs  are  superimposed  on  the  output  image  plane.  Areas  of  good  match  be¬ 
tween  the  photographs  result  in  high  fringe  contrast  in  the  output  plane. 

Since  high  metric  accuracy  is  required  in  the  contour  generating  process,  an 
array  was  chosen  to  measure  the  position  and  contrast  of  fringes  in  the  HOC 
output  plane. 

The  HOC  correlation  process  is  extremely  rapid.  The  speed  of  the  overall 
HOC  system  is  determined  chiefly  by  the  speed  at  which  image  data  can  be 
extracted  from  the  array  and  processed  by  the  computer  which  is  a  part  of  the 
HOC  system.  This  paper  describes  two  methods  for  determining  fringe  contrast 
when  a  self-scanned  array  is  used  in  the  HOC  system.  An  analysis  is  presented 
showing  the  relation  between  data  rate  and  accuracy  in  fringe  contrast  measure¬ 
ment.  Experimental  data  are  presented  showing  the  effect  of  turbulence  on 
fringe  contrast  measurement. 


The  Measurement  Problem 


Figure  1  is  a  block  diagram  of  the  HOC  system  emphasizing  the  data  acquisition 
and  processing  section.  A  Hewlett-Packard  21MX  series  computer  is  used  for 
system  control  and  data  collection  and  reduction.  From  [2]  the  light  intensity 
distribution  on  the  image  plane  due  to  the  superposition  of  the  two  aerial 
transparencies  is, 

I(x,y)  =  t^x.y)  +  t2(x,y)  +  2t1(x,y)t2(x,y)  cos  [k6(x,y)]  (1) 

where  tj  and  t2  are  the  amplitude  transmissions  of  the  transparencies,  6(x,y) 
is  the  path  length  difference  between  the  two  interferometer  arms  and 
k  =  2-rr/X, where  A  is  the  wavelength  of  the  laser  light.  If  heterodyning  is 
introduced  so  that  the  path  length  in  one  arm  is  varied  linearly  in  time  the 
intensity  becomes, 

2  2 

I  =  tx  +  t2  +  2txt2  cos  [k6  +  uit]  (2) 

where  u  is  the  angular  frequency  resulting  from  the  linear  path  length 
variation.  Each  array  element  integrates  intensity  over  its  active  area, 

A,  hence,  the  intensity  on  an  element  can  be  expressed  as, 

I  =  /tjdA  +  /t2dA  +  2/t1t2cos(k6  +  uit)dA  (3) 

The  path  length  variation  over  a  single  array  element  can  be  made  very  small 
by  adjusting  the  tilt  of  the  interferometer.  Thus,  k6  in  (3)  can  be  set 
equal  to  a  constant,  ©,  and  the  intensity  can  now  be  expressed  as, 

I  =  /tjdA  +  /t2dA  +  2cos(0  +  wt)  /tjt2dA 


(4) 


In  the  HOC  system,  image  correlation  is  determined  by  computing  a  nor¬ 


malized  correlation  coefficient  which  is  expressed  as, 
ftxt2<U< 

Cn  V/t?dA  //tadA  (5) 

This  quantity  expresses  the  degree  of  correlation  or  match  between  the  two 
transparencies  and  is  a  maximum  of  1  when  tj  =  t2  over  an  element  area.  The 
first  two  terms  in  (4)  are  the  intensities,  II  and  12  on  an  array  element  due 
to  the  photographs  in  arm  1  and  2,  respectively.  The  third  term  is  the  in¬ 
tensity  due  to  the  interference  between  the  two  photographs  and  is  directly 
proportional  to  the  correlation  between  them.  The  coefficient  in  (5)  can  be 
computed  with  the  data  from  three  measurements;  II,  12,  and  the  amplitude  of 
the  cosine  term  in  (4)  which  can  be  found  by  measuring  the  P-P  variation  of 
(4). 


The  array  used  in  the  HOC  .and  all  other  available  large  area  arrays  are 
self-scanned  devices  with  the  video  output  multiplexed  on  a  single  output 
line.  The  output  is  usually  a  train  of  pulses  or  boxcar  steps.  The  amplitude 
of  a  particular  step  is  directly  proportional  to  the  time  integrated  light 
intensity  on  the  corresponding  array  element.  The  array  is  scanned  electron¬ 
ically,  usually  in  a  raster  mode  {row  by  row),  and  the  integration  time  for 
each  element  is  the  time  required  to  scan  the  entire  array  (the  frame  period). 
Because  the  sensitivity  of  an  array  element  is  dependent  on  the  integration 
time,  self-scanned  arrays  cannot  be  randomly  accessed. 

The  measurement  of  II  and  12  with  the  array  is  straightforward,  requir¬ 
ing  only  a  single  frame  to  Input  the  steady  intensity  values  into  the 
computer.  This  can  be  accomplished  rapidly  with  a  direct  memory  access 
channel  to  the  computer.  The  measurement  of  the  peak-to-peak  value  of  (4) 


is  not  straightforward  because  the  fringe  pattern  must  be  scanned  contin¬ 
uously  across  the  array  because  of  air  turbulence.  This  situation  requires 
that  the  peak-to-peak  variation  be  found  by  sampling  continuously  the  time 
varying  pattern  and  looking  for  minimum  and  maximum  intensity  values.  The 
spatial  variation  in  the  fringe  pattern  will  be  nonuniform  in  general  so 
that  each  array  element  is  at  an  unknown  phase  point  on  the  time  varying 
sinusoid. 

The  peak-to-peak  values  can  be  determined  by  using  the  system  in  Figure  1 
in  the  following  manner.  Two  array  buffers  are  set  aside  in  the  computer 
memory;  one  for  the  minimum  values  and  one  for  the  maximum  values.  With  the 
fringe  pattern  varying  on  the  array,  one  frame  is  stored  in  both  buffers.  On 
the  next  frame  the  computer  inputs  and  compares, element-by-element, the  value 
to  the  stored  max  and  min.  If  an  input  element  value  is  less  than  the  min, 
it  replaces  the  min.  If  it  is  greater  than  the  max,  it  replaces  the  max; 
otherwise  it  is  discarded.  This  process  is  done  for  each  element, and  it  is  repeated 
for  a  predetermined  number  of  frames.  The  only  way  to  be  sure  of  the  accuracy 
of  the  max  and  min  values  is  to  take  a  sufficient  number  of  samples, yet  to  be 
determined.  The  peak-to-peak  sampling  process  is  the  bottleneck  in  the  com¬ 
putation  of  the  correlation  coefficients.  Therefore,  the  speed  at  which  map 
compilation  can  be  done  with  the  HOC  is  largely  determined  by  the  speed  at 
which  the  P-P  fringe  values  can  be  determined.  The  factors  determining  the 
speed  of  the  P-P  sampling  process  are  the  computer  speed  and  the  number  of 
frames  sampled.  The  computer  speed  can  be  maximized  by  efficient  assembly 
level  programming  or  microprogramming.  More  important  is  the  number  of  samples 
taken  to  determine  the  peak-to-peak  variations.  Here, there  is  a  trade-off  be¬ 
tween  speed  and  accuracy  and  this  sets  an  upper  bound  on  the  compilation  speed. 


The  next  section  contains  an  analysis  showing  the  relation  between  the 
number  of  samples  taken  and  accuracy  in  P-P  determination. 


P-P  Measurement  Error 


As  the  fringes  are  scanned  sinusoidally  across  the  array,  the  computer 
samples  consecutive  array  frames  to  determine  the  P-P  level  for  each  element. 
Error  arises  in  this  process  because  of  the  array  integration  time  and  be¬ 
cause  the  minimum  and  maximum  values  might  not  be  sampled  during  the  measure¬ 
ment  time.  If  Tf  is  the  array  frame  period, and  w  is  the  angular  frequency 
of  the  fringes  scanning  across  an  element,  then  the  element  integrates  the 
sinusoid  over  an  angle  a  *  wt^as  shown  in  Figure  2.  If  0  is  the  angle  to 
the  center  of  the  integration  strip  of  width  o,  the  element  response  can  be 
expressed  as 


U; 


Sine  de 


e-- 


where  £  is  a  variable  of  integration. 
When  evaluated,  (6)  becomes 


2  .  o 

—  sin  —  sine 
a  2 


or 


sine— sine 


(6) 


(7) 

(8) 


The  frame  integration  time  has  the  effect  of  reducing  the  value  measured  on 

a 

the  sine  wave  by  the  factor,  sincy.  This  in  itself  is  a  source  of  P-P  measure¬ 
ment  error.  However,  greater  error  will  occur  because  of  sampling  at  random 
points  on  the  fringe  cycle.  The  number  of  samples  per  cycle  is  N  ■  2ir/o  «  2ir /utf 
At  least  two  samples  must  be  taken  to  determine  the  P-P  value.  Since 
the  N  samples  may  be  taken  at  any  phase  on  the  fringe  cycle  the  error  will 


vary  widely.  Figure  2  shows  the  condition  resulting  in  maximum  error.  It 
is  assumed  that  an  even  integer  number  of  samples  is  taken  per  cycle.  The 
maximum  error  will  occur  when  samples  are  taken  at  equal  angles  from  90° 
but  not  at  90°.  If  the  element  integration  time  is  neglected  , 
the  maximum  error  in  percent  is. 


If  the  integration  time  is  considered  the  expression  for  maximum  error  becomes. 


r  N  *  *  1 

Em  =  100  1  -  —  sin  —  cos  — 

m  l  11  N  N  J  (10) 

Figure  3  shows  the  maximum  percent  error  in  peak-to-peak  sampling  as  a  function 
of  N,  the  number  of  samples  taken  per  fringe  cycle  (Eq.10) .  The  maximum 
error  decreases  initially  (very  rapidly)  with  increasing  N,  and  a  point  is  soon 
reached  where  very  little  is  gained  in  taking  more  samples. 


P-P  Determination  by  90°  Phase  Shift 


If  air  turbulence  is  very  small,  an  alternate  and  much  quicker  method  exists 
for  determining  the  peak-to-peak  intensity  variation  on  each  element.  If 
heterodyning  is  eliminated, the  third  term  in  equation  (4)  becomes. 


Mi  =  ^^tjtgdA 


COS  0 


(11) 


shifting  the  phase  of  one  arm  of  the  interferometer  by  90°  and  taking  another 


reading  yields. 


M2  =  [2/t1t2dA]  cos(e  +y) 

~  -  |2^txt2dAj  sine 

Thus,  the  correlation  term  in  (4)  can  be  extracted  as  below 

//v2  ■  T  J*l  *  H2 


(12) 


(13) 


Only  four  measurements  are  required  to  compute  the  normalized  correlation 
coefficient  using  this  method;  II  and  12  as  mentioned  before,  13,  (both  arms 
interfering  at  one  phase  angle),  and  14  (both  arms  interfering  at  a  90° 
phase  shift  from  13).  The  quantities  Mj  and  are  expressed  as: 


M1  =  13  -  II  -  12 
M2  =  14  -  II  -  12 


(14) 


The  phase  shift  can  be  accomplished  by  rotating  a  piezoelectric  mirror  a 
fixed  distance  between  measurements. 


Theoretically,  the  only  error  in  this  method  arises  because  of  the  PZ  dis¬ 
placement  which  can  be  made  very  small  by  calibration.  However, 
turbulence  present  in  an  actual  operating  system  may  render  this  method  im¬ 
practical  because  of  the  unknown  phase  changes  between  measurements  13  and  14. 


The  Effect  of  Air  Turbulence 

Air  turbulence  causes  random  fluctuations  in  the  fringe  pattern  when  both 


arms  of  the  interferometer  are  open.  This  random  phase  fluctuation  can  be 


modeled  as  low  frequency  noise,  N(t).  If  turbulence  is  considered,  the  third 
term  in  equation  (4)  becomes, 

2  cos[u»t  +  N(t) ]  //t1t2dA  (15) 

If  the  highest  significant  frequency  component  in  the  power  spectrum  of  N(t) 
is  small  compared  to  o>,  then  N(t)  can  be  considered  constant  over  one  cycle 
of  the  heterodyne  signal .and  the  P-P  error  analysis  holds.  Likewise,  if 
N(t)  is  essentially  constant  during  the  measurement  of  13  and  14,  the  90° 
phase  shift  method  can  be  used.  It  may  be  possible  to  shield  the  HOC  system 
so  that  turbulence  can  be  controlled  to  allow  the  use  of  the  quicker, 90°  phase- 
shift  .method.  Moderate  turbulence  can  be  tolerated  in  the 

heterodyne  method,  but  it  prohibits  the  use  of  the  90°  phase-shift  method.  Even 
if  turbulence  is  severe,  it  is  possible  to  make  measurements  with  the 
heterodyne  method  if  enough  samples  are  taken.  The  P-P  error  analysis  pre¬ 
sented  earlier  would  not  apply  in  this  case  and  the  relation  between  samples 
and  error  would  have  to  be  determined  experimentally. 


Sampling  Fringes  with  Air  Turbulence,  Experimental  Results 

Figure  3  shows  the  maximum  P-P  error  as  a  function  of  the  number  of  samples 
per  fringe  cycle  for  the  ideal  case  of  no  air  turbulence  in  the  interferometer. 

This  condition  can  be  met  only  by  shielding  the  light  paths  in  order  to  eliminate  air 
currents.  In  a  more  practical  situation  it  is  suspected  that  air  turbulence 
would  degrade  the  performance  depicted  in  Figure  3.  Figure  4  shows  the 
experimental  setup  used  to  study  the  effect  of  air  turbulence.  The  HOC  Is  used 
with  open  light  paths  to  produce  broad  fringes  on  the  array.  A  beamsplitter  is 


used  just  in  front  of  the  array- to- image  fringes  on  a  phototransistor.  This 
enables  the  frequency  of  the  fringes  on  the  Reticon  32- x  32- element  array  to 
be  determined  by  measuring  the  frequency  of  the  photo  transistor  output  with 
the  electronic  counter.  The  array-video  output  is  monitored  with  an  oscillo¬ 
scope, and  the  light  intensity  on  the  array  is  varied  with  the  neutral-den¬ 
sity  step  wedge.  The  size  of  the  fringes  on  the  array  is  monitored  with  the 
video  display.  The  video  is  digitized  with  the  A/D  converter  for  input  to  the 
computer.  The  fringes  are  scanned  across  the  array  by  varying  the  path  length 
of  one  arm  with  a  piece  of  glass  rotated  by  a  mirror  scanner.  A  triggerable 
ramp  generator  drives  the  scanner.  The  computer  is  used  for  control,  data 
reduction,  and  storage;  experimental  results  are  output  to  the  teletype. 

This  setup  was  used  to  sample  fringes  many  times  to  obtain  a  statistical 
distribution  of  sampling  error.  A  FORTRAN  program  was  written  to  take  a  spec¬ 
ified  number  of  P-P  measurements,  store  the  results, and  compute  the  probability 
of  percent  of  error.  During  the  experiment,  the  amplitude 
and  duration  of  the  ramp  was  adjusted  so  that  the  frequency  of  the  fringes 
indicated  by  the  counter ,resul ted  in  the  desired  number  of  sampl es-per-cycl e 
for  the  experiment.  For  example,  the  array  had  an  input  frame  period  of  75  ms. 
If  four  samples- per-cycle  were  desired,  the  ramp  would  be  adjusted  to  give  a 
fringe  period  of  4  x  75  =  300  ms.  Turbulence  prevented  the  exact  adjustment 
of  fringe  period.  The  ramp  was  set  to  give  the  approximate, average  period 
over  a  short  time-span.  The  helium-neon  laser  light  was  attenuated  with  the 
neutral  density  wedge  to  place  the  fringe  amplitude  within  the  dynamic  range 
of  the  array.  The  parameters  required  for  an  experimental  run  are,  (1)  the 
coordinates  of  the  array  element  where  P-P  error  is  to  be  tested,  (2)  the 
number  of  samples  per  fringe  cycle  and  (3)  the  number  of  P-P  measurement  trials. 


After  these  parameters  are  Input  with  the  TTY,  data  acquisition  and  result 
output  are  done  under  computer  control.  The  measurement  cycle  is  as  follows: 

The  computer  sends  a  trigger  signal  to  the  ramp  generator  to  initiate  fringe 
scanning.  Peak-to-peak  measurements  are  made  on  the  selected  array  element  for 
the  specified  number  of  samples  per  cycle.  The  duration  of  the  ramp  is  not 
computer  controllable,  and  a  fixed  number  of  measurements  can  be  made  before 
the  ramp  returns.  The  computer  stops  taking  measurements  before  the  end  of 
the  ramp  and  waits  for  a  period  of  time  sufficient  to  insure  the  return  of  the 
ramp  before  triggering  the  ramp  again.  This  process  is  repeated  until  the 
specified  number  of  measurements  have  been  taken  and  stored. 

The  maximum  P-P  measurement  is  taken  as  the  true  value  and  the  percent 
error  is  computed  in  terms  of  this  true  value.  The  set  of  measurements  is  divided  into 
percent  error  bins, and  the  probability  of  the  error  being  in  the  percent  range 
of  the  bin  is  computed  by  dividing  the  population  of  the  bin  by  the  total  number 
of  measurements.  Figures  5a  and  5b  show  the  probability  of  P-P  measurement 
error  for  4  and  10  samples-per-cycle  respectively.  Each  point  gives  the  prob¬ 
ability  of  the  error  being  between  the  indicated  error  and  the  next  lower  error.' 

The  maximum  theoretical  error  as  given  in  Figure  3  is  shown  in  the  plots.  The 
turbulence  results  in  a  possibility  of  errors  greater  than  the  theoretical  maxi¬ 
mum  as  the  figures  clearly  indicate.  As  expected,  turbulence  results  in  a 
greater  relative  variation  from  the  maximum  error  for  the  slower  10-sample  rate 
than  for  the  4  sample  rate.  The  effect  of  turbulence  could  be  eliminated  if  the 
computer  could  sample  at  a  sufficiently  fast  rate.  Although  air  turbulence  is 
the  major  factor,  the  figures  reflect  all  the  random  fluctuations  in  the  array 
image,  such  as  laser  intensity  variations  and  electronic  noise.  The  turbulence 
present  during  the  measurements  is  fairly  typical  of  a  laboratory  environment  con¬ 
taining  electronic  equipment  with  cooling  fans. 


Concl us Ions 


The  solid-state-sensi ng  array  is  attractive  as  a  device  for  the  measurement 
of  interferometer  fringe  contrast  and  position  because  of  its  fixed  and  accurate 
element  geometry.  Fringe  contrast  measurement  requires  the  measurement  of  peak- 
to-peak  intensity  variation  on  the  array.  If  air  turbulence  is  eliminated,  this 
can  be  accomplished  quickly  by  making  measurements  of  the  fringe  pattern  dif¬ 
fering  by  90°  phase  shift.  The  presence  of  turbulence  requires  the  random  sam¬ 
pling  of  a  time  varying  fringe  pattern  to  determine  the  peak-to-peak  intensity. 
In  this  case,  there  is  a  trade  off  in  speed  (the  number  of  samples  taken)  and 
the  maximum  allowable  error  in  peak-to-peak  measurement.  If  samples  can  be 
taken  fast  enough  to  make  the  effect  of  turbulence  insignificant,  the  maximum 
error  can  be  predicted  theoretically.  Significant  turbulence  increases  the 
likelihood  of  errors  larger  than  the  theoretical  maximum. 
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1.  HOC  Data  Processing  System 

2.  Sine  wave  sampling  at  maximum  P-P  error 

3.  Maximum  sine  wave  P-P  error  versus  number  of  samples 

4.  Experimental  Apparatus  for  P-P  sampling 

5a.  Percent  error  probability  for  four  samples  per  cycle 
with  turbulence 

5b.  Percent  P-P  error  probability  for  ten  samples  per  cycl 
with  turbulence 
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